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Abstract. We report on the fabrication of Nb/AlOrr/Pdo^Nio.is/Nb superconduc- 
tor/insulator/ferromagnetic metal/superconductor (SIFS) Josephson junctions with high critical 
current densities, large normal resistance times area products, high quality factors, and very good 
spatial uniformity. For these junctions a transition from 0- to 7r-coupling is observed for a thickness 
oIf — 6nm of the ferromagnetic Pdo.82Nio.i8 interlayer. The magnetic field dependence of the 7r-coupled 
junctions demonstrates good spatial homogeneity of the tunneling barrier and ferromagnetic interlayer. 
Magnetic characterization shows that the Pdo.82Nio.i8 has an out-of-plane anisotropy and large saturation 
magnetization, indicating negligible dead layers at the interfaces. A careful analysis of Fiske modes 
provides information on the junction quality factor and the relevant damping mechanisms up to about 
400 GHz. Whereas losses due to quasiparticle tunneling dominate at low frequencies, the damping is 
dominated by the finite surface resistance of the junction electrodes at high frequencies. High quality 
factors of up to 30 around 200 GHz have been achieved. Our analysis shows that the fabricated junctions 
are promising for applications in superconducting quantum circuits or quantum tunneling experiments. 

PACS. 74.50.+r Tunneling phenomena; Josephson effects - 74.45.+C Proximity effects; Andreev reflection; 
SN and SNS junctions - 75.30.Gw Magnetic anisotropy - 85.25.Cp Josephson devices 



1 Introduction 

The interplay between superconductivity (S) and ferro- 
magnetism (F) at S/F- interfaces and more complex S/F- 
heterostructures results in a rich variety of interesting 
phenomena [1,2,3,4,5,6]. Examples are the oscillations 
in the critical temperature of S/F-bilayers with increas- 
ing F thickness [5,7], oscillations in the critical current 
of SFS Josephson junctions [8,9,10,11,12], or the varia- 
tion of the critical temperature of FSF trilayers as a func- 
tion of the relative magnetization direction in the F lay- 
ers [13]. At an S/F-interface, the superconducting order 
parameter does not only decay in the F layer as at su- 
perconductor/normal metal (S/N) interfaces but shows a 
spatial oscillation resulting in a sign change. This oscilla- 
tory behavior is the direct consequence of the exchange 
splitting of the spin-up and spin-down subbands in the 
F layer, causing a finite momentum shift q — ±2E cx /vv 
of the spin-up and spin-down electron of a Cooper pair 
leaking into the F layer [4,5,6]. Here, E ex is the exchange 
energy and i>f the Fermi velocity in the F layer. The decay 
and oscillatory behavior of the order parameter in the F 
layer can be described by the complex coherence length 
^p 1 = £-1 +«£p 2 1 [4,5,6,14,15,16,17,18]. For the case of 
large exchange energy and negligible spin-flip scattering, 
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£fi = £f2 = \J hD I E cx with D the diffusion coefficient 
in the F layer and h = h/2w the reduced Planck con- 
stant [4,6,15]. The oscillatory behavior of the order pa- 
rameter has been unambiguously proven by the observa- 
tion of the 7r-coupled state in SFS Josephson junctions [8, 
9, 10,11,12, 19, 20, 21, 22, 23, 24, 25, 26, 27] . Here, the ground 
state has a phase difference of 7r between the macroscopic 
superconducting wave functions in the junction electrodes. 
Therefore, 7r-j unctions have an anomalous current-phase 
relation / s = I c sm(ip + tt) = —I c sin ip corresponding to a 
negative critical current —I c [3,14]. Using Josephson junc- 
tions with a step- like F layer thickness, also junctions with 
a coupling changing between and it along the junction 
have been realized [28,29,30]. 

An interesting application of 7r-coupled Josephson 
junctions are 7r-phase shifters in superconducting quan- 
tum circuits. For example, superconducting flux quan- 
tum bits [31,32,33,34,35,36] require an external flux bias 
^ext = $o(n+h) to operate them at the degeneracy point. 
Here, <Pq = h/2e is the flux quantum with the elemen- 
tary charge e and n is an integer. This requirement makes 
flux qubits susceptible to flux noise, which may be in- 
troduced through the flux biasing circuitry. Furthermore, 
to operate a cluster of (coupled) flux qubits, which in- 
evitably will have a spread in parameters, requires an 
individual and precise flux bias for each qubit. To cir- 
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cumvent this problem the insertion of 7r-phase shifters 
into the flux qubit loop has been suggested [37,38,39,40] 
and successfully demonstrated recently [41]. The appli- 
cations of 7r-coupled Josephson junctions in classical or 
quantum circuits in most cases requires high critical cur- 
rent densities J c and high critical current times normal 
resistance products, I c R n - Furthermore, they are not al- 
lowed to deteriorate the coherence properties when used 
in superconducting quantum circuits. Potential sources of 
decoherence are for instance spin-flip processes in the F 
layer or the dynamic response of the magnetic domain 
structure. Another serious difficulty comes from dissipa- 
tion due to quasiparticle excitation. It has been shown 
recently, that long decoherence times require junctions 
with both high J c and large normal resistance times area 
products, R n ■ A [42]. This is difficult to be achieved in 
SFS junctions due to the low resistivity of the metallic F 
layer. However, the situation can be improved by insert- 
ing an additional insulating barrier, resulting in a super- 
conductor/insulator/ferromagnet /superconductor (SIFS) 
stack. Here, much higher R n ■ A values at modest J c can 
be achieved. In particular, underdamped SIFS Josephson 
junctions can be realized allowing for the study of the 
junctions dynamics. Therefore, there has been strong in- 
terest in SIFS junctions recently [25,26]. 

In this article we report on the fabrication as well 
as the static and dynamic properties of a series of 
Nb/A10 I /Pd .82Nio.i 8 /Nb (SIFS) Josephson junctions 
with different thickness dp of the ferromagnetic interlayer. 
We succeeded in the fabrication of junctions with high J c 
and i? n • A values resulting in high junction quality fac- 
tors Q. We discuss the dependence of the I c R n product on 
the thickness of the F layer, including the transition from 
0- to 7r-coupling. We also address the magnetic proper- 
ties of the F interlayer in the magnetic field dependence 
of the critical current. Special focus is put on the analysis 
of Fiske modes. Comparing these resonant modes to the- 
oretical models allows us to evaluate the quality factors 
of the junctions over a wide frequency range from about 
10 to 400 GHz and to determine the dominating damping 
mechanisms. 



2 Sample Preparation and Experimental 
Techniques 

The fabrication of SIFS Josephson junctions with control- 
lable and reproducible properties was realized by the depo- 
sition of Nb/A10 x /Pd .82Nio.i 8 /Nb multilayers by UHV 
dc magnetron sputtering and subsequent patterning of 
this multilayer stack using optical lithography, a lift-off 
process, as well as reactive ion etching (PJE). Thermally 
oxidized silicon wafers (~ 50 nm oxide thickness) were 
used as substrates. 

In a first step the whole Nb/AlO^/Pdo^Nio.is/Nb 
(SIFS) multilayer stack was sputter deposited in-situ in 
a UHV dc magnetron sputtering system with a back- 
ground pressure in the low 10 -9 mbar range. This system 
is equipped with three sputter guns (Nb, Pdo.s2Nio.1s, 



Al) and an Ar ion beam gun for surface cleaning. The 
sputtering chamber is attached to a UHV cluster tool, al- 
lowing for the transfer of the samples to an AFM/STM 
system for surface characterization without breaking vac- 
uum. The multilayer stack consists of a niobium base elec- 
trode with thickness di ; Nb = 85 nm, an aluminum layer 
of thickness c?ai = 4nm, a ferromagnetic Pdo.s2Nio.1s 
(PdNi) interlayer with thickness dp ranging between 4 and 
15 nm, and finally a niobium top electrode with thickness 
<^2,Nb = 50 nm. The reproducible fabrication of ferromag- 
netic SIFS Josephson junctions requires the precise control 
of the thickness of the PdNi layer and the minimization 
of the roughness of the involved interfaces. Therefore, we 
have carefully optimized the parameters of the sputtering 
process (Ar pressure, power, substrate-target distance) for 
both the Nb and PdNi layers to obtain films with very 
smooth surfaces. After optimization, the rms roughness 
of the Nb base electrode was reduced to ~ 0.4 nm. The 
ferromagnetic PdNi layers showed a slightly larger rms 
roughness of ~ 0.8 nm. For the deposition of Nb and Al 
an Ar pressure of 2.7 x 10 -3 mbar was used, resulting in 
a deposition rate of 0.7nm/s at 200 W for niobium and 
0.2nm/s at 40 W for aluminum, respectively. PdNi was 
sputtered at a higher Ar pressure of 2 x 10 -2 mbar, result- 
ing in a growth rate of 0.4nm/s at 40 W. The transition 
temperature of the Nb films was T c = 9.2 K and the Curie 
temperature of the PdNi films was determined to about 
150 K by SQUID magnetometry. 

A critical process step in the deposition of the 
Nb/AlOz/Pdo.j^Nio.is/Nb (SIFS) multilayer stack is the 
fabrication of the tunneling barrier. To achieve good re- 
producibility, the tunneling barrier was realized by partial 
thermal oxidation of the 4 nm thick Al layer inside the 
sputtering chamber. The thickness of the AlO^ tunneling 
barrier was adjusted by varying the oxygen partial pres- 
sure and the duration of the thermal oxidation process. 
The oxidation time was varied between 60 and 240 min in 
a pure oxygen atmosphere of 0.1 mbar. We note that the 
thickness of the AlO^ tunneling barrier determines the 
R n ■ A values of the junctions, because the tunneling resis- 
tance is much larger than the resistance of the Pdo.s2Nio.i8 
layer. 

After the controlled in-situ deposition of the 
Nb/AlOz/Pdo^Nio.is/Nb multilayer stack, in the next 
step the SIFS Josephson junctions are fabricated by a suit- 
able patterning process. We used a three-stage self-aligned 
process based on optical lithography, a lift-off process and 
reactive ion etching (PJE). In the first step, the base elec- 
trode is defined by patterning a long about 100 /im wide 
strip into the whole multilayer stack. This is achieved by 
placing a suitable photoresist stencil on the Si substrate 
and using a lift-off process after the deposition of the SIFS 
multilayer stack. Fig. la shows a cross-sectional view of 
the SIFS multilayer after this step. Next, the junction area 
is patterned by etching a mesa structure into the SIFS 
stack by placing a photoresist stencil on top of the SIFS 
stack. This resist stencil serves as the etching mask in a 
RIE patterning process thereby defining the shape and size 
of the junction area. Junction areas between 2.5 x 2.5 /1m 2 



G. Wild et al.: Josephson Coupling and Fiske Dynamics in Ferromagnetic Tunnel Junctions 



3 



(a) 



(4-15nm) 



AI/AIO, 
(4 nm) 



(b) 



Nb (50 nm) 



Nb(50nm) 



Nb (85 nm) 



Si substrate 



Nb (85 nm) 



Si substrate 



(c) 



t 


Si0 2 




resist 

^ c 


\\ 


Nb (50 nm) 


[ / 








Nb (85 nm) 



(d) 



— . 



Si substrate 



□ 



Nb (85 nm) 



Si substrate 



Z3 



Fig. 1. Cross-sectional views of the SIFS Josephson junc- 
tions after different steps of the fabrication process, (a) Cross- 
sectional view after the definition of the about 100 fj,m wide 
strip forming the base electrode by a lift-off-process, (b) Cross- 
sectional view after the etching of the mesa structure defining 
the junction area by a RIE process, (c) Cross-sectional view 
after the additional oxygen plasma process reducing the size 
of the resist stencil and the subsequent deposition of the Si02 
wiring insulation. The dashed circles highlight the coverage of 
the edges of the niobium top electrode by the S1O2 layer to 
prevent shorts, (d) Cross-sectional view of the completed junc- 
tion after the deposition of the Nb wiring layer and patterning 
by a lift-off process. 



and 50 x 50 /xm 2 have been realized. The RIE process was 
performed in a SF 6 plasma (time: 70s, voltage: 300V). 
Note that the RIE process selectively patterns the Nb top 
electrode and the PdNi layer because the AlO x layer acts 
as an effective stopping layer. The resist stencil defining 
the junction area is used for the lift-off process in the 
subsequent deposition of the SiC>2 wiring insulation (self- 
aligned process). To avoid electrical shorts between the 
wiring layer for the top electrode and the base electrode, 
the lateral dimensions of the resist mask were reduced by 
about 150 nm using an oxygen plasma process in the RIE 
system immediately after the mesa patterning. In this way, 
the SiC>2 wiring insulation also covers the junction edges, 
preventing electrical shorts (cf . Fig. lc) . The 50 nm thick 
SiC>2 wiring insulation is deposited by rf-magnetron sput- 
tering in a 75%Ar/25%02 atmosphere. In a last step, the 
200 nm thick niobium wiring layer is deposited. The Nb 
deposition was done by dc magnetron sputtering and the 
patterning was realized by optical lithography and a sub- 
sequent lift-off process. To obtain a good superconducting 
contact between the Nb wiring layer and the Nb top elec- 
trode, the surface of the top electrode has been cleaned 
in-situ prior to the deposition of the wiring layer using an 
Ar ion gun. A cross-sectional view of the completed junc- 
tion is shown in Fig. Id. Of course, the junction fabrica- 
tion process allows for the fabrication of several junctions 
with different junction areas on the same wafer. In this 
way the reproducibility of the process can be checked by 
measuring the on-chip parameter spread. 



The characterization of the junctions was performed 
in a 3 He cryostat and a 3 He/ 4 He dilution refrigerator, al- 
lowing temperatures down to 300 mK and 20 mK, respec- 
tively. The dilution refrigerator was placed in a rf-shielded 
room to reduce high-frequency noise. Furthermore, exter- 
nal magnetic fields have been reduced by fj,- metal and/or 
cryoperm shields. Small magnetic fields aligned parallel 
to the junction barrier could be applied by a supercon- 
ducting Helmholtz coil. The current-voltage characteris- 
tics (IVCs) have been measured in a four-point configura- 
tion. A low-noise current source and voltage preamplifier 
(Stanford Research SR 560) have been used. Both were 
battery powered and placed inside the rf-shielded room. 
Derivatives of the IVCs have been taken using a lock-in 
technique. The magnetic properties of the ferromagnetic 
interlayer have been measured using a Quantum Design 
SQUID magnetometer. 



3 Superconducting and Magnetic Properties 

In this section we discuss the superconducting and mag- 
netic properties of the SIFS Josephson junctions fabri- 
cated with the process described above. In particular, we 
will address the magnetic properties of the Pdo.s2Nio.1s 
layer and the dependence of the product V c = I c R n of 
critical current I c and normal resistance R n of the junc- 
tions on the thickness d-p of the ferromagnetic interlayer. 
We will see that this dependence shows a clear transition 
from 0- to 7r-coupled junctions at d-p ~ 6 nm. In this arti- 
cle we will focus mainly on the properties of the 7r-coupled 
SIFS junctions. We will discuss their current-voltage char- 
acteristics (IVCs) and their magnetic field dependence of 
the critical current, allowing us to determine the basic 
parameters of the SIFS junctions. 



3.1 Magnetic Properties of the PdNi Layer 

The Curie temperature of the Pdo.s2Nio.i8 layer has been 
determined to about 150 K by measuring the magne- 
tization versus temperature dependence using a Quan- 
tum Design SQUID magnetometer. For thin ferromag- 
netic films the easy axis of the magnetization usually is 
parallel to the film plane to minimize the free energy 
contribution due to the shape anisotropy. However, in 
very thin films the magnetic surface energy may be dom- 
inant resulting in an easy axis perpendicular to the film 
plane. To get information on the direction of the easy axis 
we have recorded the magnetization versus applied mag- 
netic field curves of the thin Pdo.82Nio.i8 layer within a 
Nb/A10a;/Pdo.82Nio.i8/Nb multilayer stack with the mag- 
netic field applied in- and out-of-plane. The result is shown 
in Fig. 2. Qualitatively, an almost rectangularly shaped 
M(H) curve is expected for the field applied along the 
easy axis because the magnetization tries to stay along this 
preferred direction as long as possible and than abruptly 
switches to the opposite direction at the coercive field. 
In contrast, for the field applied along the hard axis the 
magnetization is expected to be gradually rotated out of 
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Fig. 2. Magnetization versus applied magnetic field 
curves of a 8.4 nm thick Pdo.82Nio.i8 layer sandwiched in a 
Nb/AlOz/PdNi/Nb multilayer stack for the magnetic field ap- 
plied parallel (top panel) and perpendicular (bottom panel) to 
the film plane. The curves were measured at 11 K, that is, just 
above the critical temperature of the Nb films. The magnetiza- 
tion data have been corrected by subtracting the diamagnetic 
background due to the metallic layers and the substrate. This 
contribution is determined from the slope of the M (H) curves 
at high fields above 1 T where the magnetization of the ferro- 
magnetic layer is fully saturated. 



the easy axis direction into the hard axis direction on in- 
creasing the applied magnetic field, resulting in a gradu- 
ally increasing and decreasing magnetization when sweep- 
ing the field. As demonstrated by Fig. 2, the measured 
M{H) curve is much more rounded and has smaller re- 
manent magnetization, when the field is applied in-plane, 
whereas it has an almost rectangular shape when the field 
is applied out-of-plane. This is clear evidence for an out- 
of-plane anisotropy of the PdNi film. That is, the out- 
of-plane (in-plane) direction is the magnetic easy (hard) 
axis in agreement with literature [27]. An out-of-plane 
anisotropy has been observed also for CuNi [43]. A more 
detailed analysis of the magnetic anisotropy of the PdNi 
films would require additional experiments such as ferro- 
magnetic resonance and methods providing information 
on the domain structure. However, it is difficult to per- 
form such experiments with ferromagnetic layers enclosed 
by two metallic Nb layers. 

The saturation magnetization M s of the Pdo.s2Nio.1s 
layer is about /zoM s ~ 0.2 T in good agreement with lit- 
erature values [27]. This magnetization corresponds to a 
magnetic moment of slightly above 1 /ib per Ni atom, if 
one assumes that there is negligible contribution of Pd. 
Here, is Bohr's magneton. This is close to the values 
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Fig. 3. Dependence of the / c i? n -product of SIFS junctions 
on the thickness <1f of the ferromagnetic Pdo.82Nio.i8 layer 
for samples with different junction areas. The broken lines 
are guides to the eye. The olive line is a fit of the data 
to the theory of Buzdin et al. [15] using £f = 3.88 nm and 
TTA 1 A 2 /ek B T c = 280 ^V. 



in polycrystalline bulk samples, for which /iNi — 1-1 Mb 
and /ipd — 0.1 /^B has been reported [44]. This suggests 
that there are no significant magnetic dead layers. Since 
M g was measured with an only 8 nm thick PdNi film, even 
very thin dead layers at the interfaces would result in a 
significant reduction of M s . In the measurement of the 
magnetic field dependence of the critical current of the 
SIFS junctions only an in-plane magnetic field of the or- 
der of a few mT is applied. In this small field range only 
the virgin M(H) curve is relevant, which is about linear 
and almost non-hysteretic. From the experimental data 
taken at 11 K the slope of the M(H) curve is determined 
to m = dM/dH = 1.8. 



3.2 Crossover between 0- and 7r-Coupling 

We next discuss the dependence of the characteristic junc- 
tion voltage V c = I c R n on the thickness dp of the ferro- 
magnetic interlayer. The V c values of a large number of 
junctions with different dp and junction areas A ranging 
from 5x5 /xm 2 to 50 x 50 /Ltm 2 are plotted in Fig. 3 ver- 
sus the thickness d? of the F layer. To achieve this data, 
a series of junctions with different dp was fabricated with 
otherwise identical parameters. In particular, all junctions 
have AlO x barriers obtained by a 90 min long thermal ox- 
idation process, resulting in R D ■ A ~ 40 i?/im 2 . This series 
clearly shows a change of sign in the slope of the I c R n (dp) 
dependence at a Pdo.s2Nio.1s layer thickness of g?f ~ 6 nm. 
At this dp value I c approaches zero. This feature is a clear 
signature of the crossover from 0- to 7r-coupling on increas- 
ing dp, corresponding to a change of sign of I c . Obviously, 
in the experiment we can only measure the modulus of 7 C . 

In order to theoretically describe the behavior shown 
in Fig. 3 one has to distinguish different regimes defined 
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by three energy scales [4]. These scales are the exchange 
energy E cx in the ferromagnet, the energy gap A of the 
superconductor, and h/r, where r is the elastic scatter- 
ing time in the ferromagnet. In our samples, E cx ^> A 
(A — 1.5 meV for Nb) as in the overwhelming part of 
other reports on SFS or SIFS junctions. However, there is 
a significant variation in H/t relative to the other two en- 
ergy scales. The true clean-limit holds for E cx , A ^> H/t. 
In this case the mean free path £ — vpr is large com- 
pared to the clean-limit superconducting coherence length 
£ s = hvp/irA and exchange length £p = fivp/2E cx . Here, 
vp is the Fermi velocity in the respective material. The 
true dirty limit holds for E cx , A <C H/t. In this case the 
mean free path I = vpr is the smallest length scale and the 

dirty-limit superconducting coherence length £ s c 



^KD/ttA and exchange length £p — \J £,f£ — \/foD/E cx 
are given by the geometric mean. Here, D is the diffu- 
sion coefficient in the respective material. There is also 
an intermediate regime, where E ex > H/t and A < H/t, 
which is the most complicated situation. This regime may 
particularly apply for a ferromagnet with large E ex . 

Due to the small mean free path in the PdNi alloy and 
Nb films, for our samples the simple dirty limit holds. 
Several theoretical models have been proposed for this 
limit [4, 9, 14, 15, 16, 17, 18] . For weak ferromagnets such as 
PdNi the spin-up and spin-down subbands can be treated 
identically (same Fermi velocity and mean free path) re- 
sulting in a single characteristic length scale 



hD 

E cx 



(1) 



for the decay and oscillation of the critical current as a 
function of dp. We note, however, that in the presence of 
spin-flip or spin-orbit scattering the decay and oscillation 
of I c is governed by two different length scales [14, 16, 17, 
18]. 

The solid olive line in Fig. 3 is obtained by fitting the 
data using the simple dirty limit expression of Buzdin et 
al. [15]. The fit yields £ F = 3.88 nm and = 280 pV. 

Here, T c — 9.2 K is the critical temperature of the nio- 
bium layers and A\ and Ai are the superconducting order 
parameters just at the boundary with the ferromagnetic 
layer. They are certainly much smaller than the super- 
conducting gap in bulk niobium but difficult to be de- 
termined in the geometry of our experiment. The value 
£f = 3.88 nm found for our SIFS junctions agrees very 
well with literature values [9,15,45]. The value of 280 /iV 



ttA 1 A 2 
ek B T c 



for our SIFS junctions is larger than the 



found for 

value of 110 /iV obtained by Kontos and coworkers [9,15]. 
Unfortunately, the detailed comparison of different exper- 
iments is difficult because the authors often do not state 
whether they are plotting I c R n (dp) or I c R ss (dp). If we use 
the much higher subgap resistance i? sg , the corresponding 
value derived for ""f 1 i? 2 would be about ten times larger. 
We finally note that the measured I c (dp) dependence can 
be well explained by dirty limit theory with a single length 



scale £p, suggesting that spin-flip or spin-orbit scattering 
do not play a dominant role. 

Beyond the parameter £p describing the characteristic 
length of superconducting correlations in the F layer, the 
dimensionless parameter 



7b = 



PB&F 

IT 



(2) 



is used to describe ferromagnetic Josephson junctions. It 
characterizes the transparency of the F/S interfaces with 
Pb the interface resistance times area and ctf the conduc- 
tivity of the F layer. In general, two 7b values for the 
two F/S interfaces have to be used. In our experiment 
the presence of the additional AlO^ barrier at one S/F 
boundary can be modeled by a very low transparency in- 
terface (7b i ^> 1), while the other boundary has high 
transparency (7B2 <C 1) for the in-situ fabricated stacks, 
i.e. 7b — 7bi- We further note that in general the mea- 
sured total R n ■ A product can be expressed as 



R n ■ A = p tu 



Pint 



dp 



(3) 



Here, p t un and p mt are the resistance times area val- 
ues due to the tunneling barrier and the PdNi/Nb in- 
terfaces. For our junctions, ptun 3> Pint and, moreover, 
Ptun ^ dp/ up. With ct pdNi — 10 7 i7 _1 m _1 we estimate 
dp/ap ~ 10 _3 i?/im 2 which is by about five orders of mag- 
nitude smaller than the measured R n ■ A values. That is, 
the contribution of the F layer to R n ■ A is negligible. 
Hence, R n ■ A ~ pt un , meaning that the measured R n ■ A 
values are dominated by the AlO^ tunneling barrier as ex- 
pected. In this case we can write 7b — R n Aap/^p. With 
CTpdNi — 10 7 i7 _1 m _1 , £p = 3.88 nm and the measured 
R n ■ A value of about 40 Qpm 2 for this junction series we 
estimate 7b — 10 5 . This high value is not surprising due 
to the additional tunneling barrier in our junctions. 

The derived £f values can be used to estimate the 
exchange energy in PdNi. Using a Fermi velocity vp ~ 
5 x 10 5 m/s [46] and the fact that the mean free path £ 
for the very thin PdNi layers is about given by the film 
thickness dp, we obtain D ~ vp£ ~ up^F = 5 x 10 _3 m 2 /s 
for dp — 10 nm. With this value we derive E cx ~ hD/£ F ~ 
20mcV using £p = 3.88 nm. This value is in good agree- 
ment with values between about 10 and 50 meV quoted 
in previous work [9,21,27,47]. With the same numbers we 
obtain H/t ~ hvp/dp ~ 30meV. This shows that the SIFS 
junctions with PdNi interlayer are close to the intermedi- 
ate regime since E cx ~ H/t. 



3.3 IVCs and Magnetic Field Dependence of the 
Critical Current 

In the following we discuss the IVCs and magnetic field 
dependence of the critical current. We focus on SIFS junc- 
tions with a F layer thickness of dp = 8.4 nm resulting in 
7r-coupling. The AlO^ tunneling barrier has been achieved 
by thermal oxidation for 4 h in pure oxygen to obtain high 
R n ■ A values and, in turn, low damping of the self excited 
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resonances. We start our discussion by considering the ef- 
fect of the specific junction geometry (cf. Fig. 4) used in 
our experiments on the derived junction parameters. Al- 
though these effects are often neglected in literature, it is 
mandatory to take them into account in a detailed evalu- 
ation. Of particular importance are the finite thickness of 
the junction electrodes and the so-called idle region. The 
latter is formed when the wiring of the top electrode is de- 
posited to complete the junction structure. Then the bot- 
tom electrode, the SiC>2 wiring insulation and the wiring 
layer are forming an SIS structure next to the junction 
area. 



Effective Magnetic Thickness It is well known that a 
magnetic field applied parallel to the surface of a bulk 
superconductor decays exponentially inside the supercon- 
ductor due to the Meifiner effect [48]. The characteris- 
tic screening length is the London penetration depth Al 
which is about 90 nm for Nb [49]. However, in our junc- 
tions we are using electrodes of finite thickness d. In such 
a thin film superconductor the magnetic field dependence 
perpendicular to the film (z-direction) is obtained to 

rrf \ _ #cxt,l + H cyit .2 CQSh(,Z / 'A L ) 

[Z) ~ 2 cosh(d/2A L ) 

H cxt ,i ~ H cxt:2 sinh(z/A L ) 

2 sinh(d/2A L ) ' 1 ' 

by solving the London equations [48]. Here, the film is 
assumed to extend from — d/2 to +d/2, and H cxt ,i and 
H C xt,2 are the external magnetic fields applied parallel 
to the film at both sides. The boundary conditions and 
the resulting current distributions are different for dif- 
ferent physical situations. If the field is confined on one 
side of the superconducting film, we have H cxt ,i = Ho 
and -ffoxt.2 = 0. For Josephson junctions this situation ap- 
plies to cases where one is considering only the fields re- 
lated to supercurrents flowing in the junction electrodes. 
If d < Al, screening currents are confined to a length 
scale smaller than Al, resulting in an enhanced kinetic 
inductance L s = /ioAl coth(d/AL) as compared to a bulk 
superconductor with L a = [IqXl [50]. That is, the thin 
film superconductor behaves equivalent to the bulk one 
with an effective screening length Al coth(rf/AL). Then, 
the magnetic penetration in the barrier layer and the junc- 
tion electrodes can be described by the effective magnetic 
thickness 

t J B = f + fidp + Ali coth + Al2 coth . (5) 

Ali Al2 

Here, t s is the thickness of the oxide barrier, dp the thick- 
ness of the F layer, e?i and c?2 are the thicknesses of the 
Nb junction electrodes, and Ali and Al2 the corresponding 
penetration depths. For our junctions we have t 3 = 4nm 
(we assume for simplicity that the oxide thickness is equal 
to the thickness of the Al layer) , d\ = 85 nm, d 2 — 250 nm, 
Ali = Al2 = ANb = 90 nm [49], g?f = 8.4 nm and fi = 1.8 
resulting in t& — 230 nm. For the idle region next to the 




Fig. 4. Sketch of the SIFS junction geometry. The junction 
area of length U and width W 1 is formed by the junction elec- 
trodes of thickness di and d^ separated by the tunneling barrier 
with thickness t 1 and dielectric constant e J and the ferromag- 
netic layer with thickness dp. Next to the junction area there 
is an overlap of the Nb wiring layer and the base electrode, 
forming the so-called idle region of length Z//2 at each side 
(1/ = a few 10 ^m and W 1 ~ in our junctions). Both layers 
are separated by the SiC>2 wiring insulation of thickness t 1 with 
dielectric constant e 1 . 



junction area (cf. Fig. 4) an expression equivalent to eq.(5) 
is obtained by replacing V by t 1 on the right hand side and 
setting dp = 0. 

Another interesting case is the situation where the 
same magnetic field is applied at both sides of the su- 
perconducting film. Here, H cxt ,i — H cxt ,2 = Hq and 
hence H(z) = H cos1i(z/Al)/ cosh(d/2AL). This case ap- 
plies when we are measuring the dependence of the junc- 
tion critical current on a magnetic field applied paral- 
lel to the junction electrodes. Since the junction (we are 
considering only short junctions here) cannot screen the 
applied field from the region between the junction elec- 
trodes, the same field is present at both surfaces of the 
junction electrodes. To derive an effective magnetic thick- 
ness for this situation we consider the total flux thread- 
ing the junction. The latter is obtained by integrating 
eq.(4) along the z-direction. Whereas J H(z)dz = 2iJ AL 
for bulk electrodes (d 3> Al), for thin film electrodes we 
obtain J H(z)dz = 2H a \\^ tanh(<f/2AL). That is, regard- 
ing the flux the thin film superconductor behaves equiv- 
alent to the bulk one with an effective screening length 
Al tanh(c?/2AL). Then, for the description of the total flux 
threading the junction we can use the effective magnetic 
thickness [51] 

? B — ^ + ^d F + A L i tanh — h A L 2 tanh -y 2 - . (6) 

2Ali 2Al2 



Effective Swihart Velocity and Frequencies of Resonant 
Modes The idle region shown in Fig. 4 influences the dy- 
namic properties of the junction. The reason is that the 
idle region acts as a dispersive transmission line parallel 
to the junction edges. As discussed in detail in section 4, 



G. Wild et al.: Josephson Coupling and Fiske Dynamics in Ferromagnetic Tunnel Junctions 



7 



the junction forms a transmission line resonator support- 
ing so-called Fiske resonances [52]. For a junction with- 
out any idle region these resonances appear at frequencies 
{jJnjl-K = n{c/2D) corresponding to voltages 



V n = n 



2D 



Here n = 1, 2, 3, . . . and 



e 3 ti 



(7) 



(8) 



is the Swihart velocity [50] with c the velocity of light in 
vacuum. Evidently, c<cas the junction barrier thickness 
t 3 is by almost two orders of magnitude smaller than the 
effective magnetic thickness i J B of the junction. 

For a junction with idle region, the resonance frequen- 
cies and corresponding voltages are shifted. We first dis- 
cuss the effect of a lateral idle region extending paral- 
lel to the resonant mode. In this case we have to con- 
sider the junction transmission line in parallel with the 
two transmission lines of the idle regions at both sides of 
the junction. To each isolated transmission line we can 
assign the inductance C 1 ' 3 — fiol^/W 1 ^ and capacitance 
C 1 ' 3 = eoe^W 1 ' 3 /t 1 ' 3 each per unit length, where the indices 
i and j refer to the idle and junction region. The corre- 
sponding phase velocities are = v C^C 1 ' 3 with v ph 



c. 

The phase velocity of the combined transmission lines is 
given by v ph = 1/y/CC with Cr 1 = (Oy 1 + {C 3 )' 1 and 
C = C 1 + C 3 . Evidently, the idle region increases the ca- 
pacitance and decreases the inductance per unit length. 
However, the inductance effect is dominant since usually 
*bAb t 3 /t 1 - This is true also for our samples. Therefore, 
the phase velocity is increased by the idle region. A more 
detailed analysis yields [53,54,55] 



Vph 



\ 



1 



Hi 



1 



(9) 



We next discuss the effect of a longitudinal idle region 
extending perpendicular to the resonant mode. In this sit- 
uation the idle region can be considered as a lumped ca- 
pacitance loading the junction at its ends, since the wave 
length usually is much larger than the dimension of the 
idle region. For short junctions (L 3 ,W 3 <C Aj) detailed 
calculations yield [55] 



1 



<J p h 



1 



C 



1-1 



(10) 



Here, C\ is the lumped load capacitance and Co the total 
junction capacitance. We see that the effect of C\ is to 
decrease the phase velocity. However, since usually t'/e 1 ^> 
t 3 /e 3 we have C\ <C Co, that is, the effect is quite small 
even if the idle region has similar area as the junction. 
With i 1 = 50 nm, e ! ~ 3.9, t 3 = 2nm, e 3 ~ 9.1, t 3 B ~ 
230 nm, ~ 260 nm and W l /W 3 = 0.8 (typical for a 
50 x 50 /jm 2 junction) we estimate v p h ~ 1.22 c for the 
lateral and u p h ~ 0.95 c for the longitudinal mode. 



(a) 



-2 



"200" 



100 



T=40mK 




A= 20x20 \irrf 
T=2K 




/ 120 



/ < 80 
/ 3, 



- 40 







\ T=500mK 


: / 






Fig. 5. IVCs obtained for 7r-coupled SIFS junctions with 
A = 20 x 20 /im 2 (a) and 50 x 50 /im 2 (b) . The dashed lines indi- 
cate the ohmic behavior V = IR n approached at large voltages. 
The upper and lower insets show enlarged views of the IVCs 
around zero voltage and the magnetic field dependence of I c , 
respectively. The magnetic field is applied parallel to the junc- 
tion barrier. The total magnetic flux $ — $ h + threading 
the junctions originates from the applied magnetic field ($h) 
and the magnetization of the F interlayer ($m). The IVCs are 
measured at 2K, while the / c (^) dependencies are taken at 
500 mK. The red IVCs in the top inset of (a) is measured at 
40 mK. 



Current- Voltage Characteristics Figure 5 shows the 
IVCs of two 7r-coupled SIFS junctions with junction ar- 
eas A = 20 x 20 fim 2 and 50 x 50 ^m 2 . From the IVCs 
we can determine several relevant junction parameters. 
First, the current value for the switching from the zero 
to the finite voltage state and vice versa gives the criti- 
cal current I c and the so-called retrapping current I r (see 
insets of Fig. 5). The corresponding current densities J c 
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and J x are obtained by dividing by the junction area A. 
From J c together with the effective magnetic thickness 
t B of the junction (cf. eq.(5)) the Josephson penetration 



depth Aj = ^h/2e/j, t 3 B Jc is derived. At T = 2K we ob- 
tain Aj ~ 70 /im. Since this value is larger than the lateral 
junction dimensions, we are dealing with small Joseph- 
son junctions. Using the simple resistively and capacitively 
shunted junction (RCSJ) model [56], we can further de- 
rive the junction quality factor Qivc — — W p/ W RC = 
4J c /7r/ r . Here, /3c is the Stewart-McCumbcr parameter, 
uj p the junction plasma frequency, and 1/wrc the RC 
time constant of the junction. Second, from the asymp- 
totic behavior at voltages large compared to the gap sum 
voltage the normal resistance R D and the normal resis- 
tance times area product, R n ■ A, are obtained which are 
about temperature independent. The R n value has to be 
distinguished from the temperature dependent subgap re- 
sistance R sg obtained from the slope of the IVCs at volt- 
ages well below the gap sum voltage. For SIFS junctions, 
i? sg is expected to increase with decreasing temperature 
in agreement with the experimental data. At 2K, for our 
SIFS junctions R sg is almost an order of magnitude larger 
than R n . 

We note that the values of I c and I T may be reduced 
and enhanced, respectively, by premature switching due to 
thermal activation or external high-frequency noise. This 
is particularly true for small area junctions with small ab- 
solute values of I c and I r . In turn, this results in reduced 
values of the quality factor. This effect is shown in the 
inset of Fig. 5a, where the IVCs of a 20 x 20 /mi 2 junc- 
tion are shown for T = 2K and 40 mK. The 40 mK data 
were taken in a well shielded dilution refrigerator using 
various filters at different temperature stages in the cur- 
rent and voltage lines, including stainless steel powder fil- 
ters. Clearly, at 40 mK significantly larger I c and smaller 
J r values are observed resulting in an about three times 
larger quality factor Qrvc = 6.2. These enhanced/reduced 
values are only partly caused by lowering the temperature 
but mostly by the reduced thermal and external noise. For 
large area junctions this effect is negligible. Here, J c and I r 
are large, making the relative effect of the equivalent noise 
current very small. For example, at T = 2K Qrvc = 5.7 is 
obtained for the 50 x 50 /im 2 junction, whereas Qivc < 2 
for the 20 x 20 /im 2 junction fabricated on the same chip. 
This small quality factor is related to the smaller I c and 
J r values of the small area junction, making it more sus- 
ceptible to thermal and external noise. In Table 1 we have 
tabulated the junction parameters derived from the IVCs 
of the 50 x 50 /im 2 and 20 x 20 /im 2 junction at 2 K and 
40 mK, respectively. 



Magnetic Field Dependence of the Critical Current The 

bottom insets in Fig. 5 show the magnetic field depen- 
dence of the critical current of the SIFS Josephson junc- 
tions measured at 500 mK. The dependencies are close to 
a Fraunhofer diffraction pattern 



Table 1. Critical current and resistance values obtained from 
the IVCs of the SIFS junctions with dp = 8.4 nm shown in 
Fig. 5 and parameters derived from them. The quantities are 
listed for a 50 x 50 /im 2 and 20 x 20 /im 2 junction at 2K and 
40 mK, respectively. 



A (/im 2 ) 


20 x 20 


50 x 50 


T (K) 


40 mK 


2K 


ic (M) 


131 


555 


Jc (A/cm 2 ) 


33 


22 


Ir (fJ,A) 


27 


125 


Rn (O) 


0.33 


0.051 


R n ■ A (Ofim 2 ) 


133 


128 


R ss (fl) 


2.14 


0.44 


I C Rn (jlV) 


44 


28.3 


I c R sg (jjY) 


280 


244 


Aj(/xm) 


59 


71 


Qivc 


6.2 


5.7 



expected for an ideal short Josephson junction with a spa- 
tially uniform J c . This demonstrates that our SIFS junc- 
tions have good uniformity of J c across the junction area, 
that is, a spatially homogeneous tunneling barrier and fer- 
romagnetic interlayer. In particular, there are no short cir- 
cuits at the junction edges or the surrounding Si02 wiring 
insulation. We note that direct information on junction 
inhomogeneities on smaller length scales can be obtained 
by Low Temperature Scanning Electron Microscopy [57, 
58,59,60]. 

We note that the total flux <P threading the junction 
is composed of the flux # H due to the applied external 
magnetic field and the flux <Pm due to the magnetization 
of the ferromagnetic layer. The two components are given 

by 



<P B = t B V fi H , 
<P M = d F V fx M . 



(12) 
(13) 



Here, U is the lateral dimension of the rectangular junc- 
tion perpendicular to the field direction, dp the thickness 
of the F interlayer, and r B the effective magnetic thickness 
of the junction given by eq.(6). Since the applied mag- 
netic field H$ = <t>a/ fAoUtr-g required for the generation 
of a single flux quantum in the junction is less than about 
1 mT for 10 /im < I) < 50 /im, the typical field range used 
in the measurement of the I c (^/<?o) curves of Fig. 5 is re- 
stricted to less than about 10 mT. For such small in-plane 
magnetic fields the M(H) curve of the ferromagnetic PdNi 
layer (cf. Fig. 2) can be well approximated by a linear de- 
pendence M ~ /j,H with /i = 1.8. With this approximation 
the total magnetic flux threading the junction can be ex- 
pressed as 



I C (H) = / c (0) 



sin(7T^/^ ) 



ir<P/<P 



(11) 



$ = H H L>(P B -d F ) 



= noH L% 1 + ( M - 1) 




(14) 
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Here, M is the magnetization component parallel to the 
applied magnetic field. If there is a significant in-plane 
magnetic anisotropy of the ferromagnetic material, this 
component may be much smaller than the absolute value 
of the magnetization. Furthermore, there may be a com- 
plicated domain structure. In this case, M is the average 
magnetization parallel to the field direction. 

We can use (14) to estimate the magnitude of the ad- 
ditional flux due to the F layer. For the junctions 
of Fig. 5, dp = 8.4 nm and hence dp/i^ ~ 0.05. There- 
fore, for fi — 1.8 the second summand in the brackets 
of (14) amounts to only about 0.04. That is, compared 
to a junction without F interlayer the total flux is en- 
hanced only by about 4%. Due to the uncertainties in 
Al and the geometrical dimensions of the junctions, this 
small effect is difficult to prove. Furthermore, since the 
virgin part of the hard axis M{H) curve of the PdNi layer 
is about linear at small fields with negligible hysteresis, 
the 2c(<£/^o) curves are expected to show negligible hys- 
teresis on sweeping back and forth the applied magnetic 
field. This is in agreement with our data and measure- 
ments on SIFS junctions with ferromagnetic CuNi inter- 
layers [61]. We note, however, that <?m strongly depends 
on magnetic history. After having increased the applied 
field to large values saturating the F layer, there is a signif- 
icant remanent magnetization on reducing the field to zero 
again. This remanence causes a large which may ex- 
ceed 100<?o- Unfortunately, applying high magnetic fields 
also results in the trapping of Abrikosov vortices in the 
Nb junction electrodes, making the interpretation of the 
measured J c (#/#o) dependencies difficult. 



4 Fiske Resonances 

According to the Josephson equations, the supercurrent 
across a Josephson junction oscillates at a constant fre- 
quency wj = 2ttV/ < Pq when a constant voltage V is applied 
across the junction [62]. On the other hand, the junction 
geometry forms a transmission line resonator of length D 
(cf. Fig. 4) with resonance frequencies uJn/2-K — rc-ofy- 
Here, n = 1, 2, 3, . . . and c is the Swihart velocity [50]. At 
an applied junction voltage V n — u} n <P /27r — n^jf, the 
oscillation frequency of the Josephson current matches the 
n-th harmonic of the junction cavity mode, potentially 
resulting in the excitation of the cavity modes. The ex- 
citation of the cavity resonances is most effective, if the 
spatial period of the Josephson current distribution along 
the junction is about matching the spatial period of the 
n-th resonant electromagnetic mode of the junction. Since 
for short junctions the Josephson current density is spa- 
tially uniform at zero magnetic field, there is no excita- 
tion of the resonant modes. However, a spatial modula- 
tion can be easily achieved by applying a magnetic field. 
The resonance is a highly nonlinear process which results 
in self-induced current steps at V n in the IVCs denoted 
as Fiske steps [52,63,64]. These steps are clearly visible 
in the IVCs shown in Fig. 6 and their derivatives dV/dl 
plotted in Fig. 7. The data were obtained for a junction 




V(uV) 

Fig. 6. Current-voltage characteristics (IVCs) of a SIFS 
Josephson junction with A = 50 x 50 fim 2 and dp = 8.4 nm 
measured at 500 mK for different magnetic fields applied paral- 
lel (bottom) and perpendicular (top) to the bottom electrode. 
The first three Fiske resonances are clearly seen and labeled 
Fl, F2, and F3. 

with area A = 50 x 50 /im 2 recorded at different mag- 
netic fields applied parallel (x-direction) and perpendic- 
ular (y-direction) to the bottom electrode. The dV/dl 
vs. V curves were measured by a lock-in technique. In 
the case of negligible damping the resonances are very 
sharp at almost constant voltages. However, in reality the 
resonances are damped due to different loss mechanisms. 
For planar type SIFS junctions as studied in our work 
the junction damping due to radiation losses is small due 
to the large electromagnetic impedance mismatch at the 
junction boundaries. Then, the quality factors Q n of the 
Fiske resonances are given mainly by the internal losses, 
most likely due to quasiparticle tunneling and the finite 
surface resistance. In general, the detailed analysis of the 
Fiske resonances can provide information on the damping 
mechanisms in Josephson junctions at very high frequen- 
cies. 

When analyzing the voltage position of the Fiske res- 
onances in detail we have to take into account the idle 
region next to the junction area. As discussed above, this 
idle region results in an increase or decrease of the phase 
velocity v p h compared to the Swihart velocity c of an ideal 
junction without any idle region, depending on whether 
the idle region is lateral or longitudinal, respectively. In 
turn, this causes an increase or decrease of the character- 
istic voltages V n = n -fffi ■ For H\\y, the resonant modes 
are extending in rr-direction, that is, parallel to the idle 
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region (lateral mode). In this case the effect on the phase 
velocity is described by eq(9) and we expect a slight in- 
crease of w p h and V n . Analogously, for H\\x the resonant 
modes are extending in y-direction, that is, perpendicu- 
lar to the idle region (longitudinal mode). In this case, 
the effect on the phase velocity is described by eq(10) 
and we expect a slight decrease of w p h and V n . This is 
in good qualitative agreement with our experimental ob- 
servation. For H\\x, resonant modes are found at 167 /iV, 
322 ^V, 475 [N, 630 fiV and 772 /iV, whereas for H\\y the 
resonant modes appear at slightly larger voltages 181 /iV, 
336 fiY, 483 ^V, 636 /iV and 780 //V. We note, that the 
correction factors estimated from eqs.(9) and (10) indi- 
cated an even bigger difference between the lateral and 
longitudinal mode. However, the quantitative evaluation 
of this difference depends on the details of the junction 
geometry and properties of the involved materials (e.g. 
dielectric constants) and therefore needs a more elaborate 
effort [53,54]. Qualitatively we can say that the effect of 
the idle region is reduced on going to higher frequencies 
(smaller wavelengths) due to a stronger confinement of the 
modes in the junction area. That is, the difference of the 
resonance voltages of the lateral and longitudinal mode 
is decreasing on going to higher harmonics. This nicely 
agrees with our observation and reports in literature [53]. 

The Fiske steps AI n (V, •P) are obtained from the mea- 
sured IVCs (cf. Fig. 6) by subtracting an ohmic back- 
ground which is determined by the about constant subgap 
resistance. Equivalently, AI n (V, #) can be obtained by in- 
tegration of the dV I dl curves shown in Fig. 7, again after 
substraction of the subgap resistance. The latter method 
yields better results for the higher modes. A detailed quan- 
titative description of the voltage position and the shape of 
the Fiske resonances was given by Kulik for standard SIS 
Josephson tunnel junctions [65] . In this section we analyze 
the Fiske steps AI n (V, <P) in the IVCs of our SIFS Joseph- 
son junctions using an extension of Kulik's theory [64,66, 
67,68]. According to this theory, the Fiske steps can be 
expressed as [64,66] 



AI n (V, 



4ir 2 7 



E 



l/<3, 



*=°Hf) 2 ] +(&) 



F 2 M) ■ (15) 



Here, Q n is the quality factor of the n-th resonant mode 
and V n its voltage position. The function F n ((p) describes 
the flux dependence of the resonant mode given by 



F„(4>) 



20COS(7T</)) 

7T [(t> 2 - (n/2) 2 ] 

20sin(7r0) 
7T [</) 2 - (n/2) 2 } 



for n = 1, 3, 5, . . . (16) 
for n= 2,4,6,... (17) 



with <p = <P/<Po. The flux dependence of the maximum 
height of the Fiske steps is given by [64,66] 



AL 



max / : max 



) = h 



Qn 

4:Tr 2 ri 




F 



c )- (18) 



Fig. 7. Differential resistance dV/dl vs. applied voltage V 
measured with a lock-in technique at 500 mK for different 
magnetic fields applied parallel (bottom) and perpendicular 
(top) to the bottom electrode. The data is taken simultane- 
ously to those shown in Fig. 6. Fiske resonances are observed 
up to almost 800 /iV corresponding to a frequency of about 
400 GHz. The origin of the resonance labeled HF at half the 
first Fiske voltage is discussed in the text. The thick red lines 
mark the measurements taken at zero applied magnetic field. 
These curves are not shifted vertically, whereas all other curves 
are subsequently shifted by 0.1 fi for clarity. 



Fitting the current steps AI n (V) measured at fixed ap- 
plied flux to (15) allows us to determine the voltage po- 
sitions V n _\\ (V nt ±) and the maximum step heights AI™^ X 
(AI™ la £) for the magnetic field applied parallel (perpendic- 
ular) to the base electrode. For low n, the quality factors 
determined by eq.(18) are not within the validity of Ku- 
lik's theory applicable only for low Q. Therefore, the the 
quality factors Q„.|| {Q n ,±) listed in Table 2 have been cal- 
culated by an extended theory more appropriate for the 
high-Q regime [64]. 

As discussed above, for short Josephson junctions the 
excitation of the cavity resonances is most effective, if the 
spatial period of the Josephson current distribution along 
the junction is about matching the spatial period of the n- 
th resonant electromagnetic mode of the junction. There- 
fore, the height of the Fiske steps strongly depends on 
the applied magnetic field and has a pronounced maxi- 
mum at a particular field value. This is shown in Fig. 8, 
where we have plotted the height Z\/™ ax of the n-th Fiske 
step together with the critical current J c of a 50 x 50 jxm 2 
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Table 2. Voltage position V n ,\\ (V n ,±) and step height AI™* (AI™X) °f tne n_tn Fiske steps obtained for the magnetic 
field applied parallel (perpendicular) to the bottom electrode. The data was derived from the IVCs of a 50 x 50 /im 2 junction 
measured at 500 mK. Also listed are the quality factors Q n ,\\ and Qn,± of the n-th Fiske resonance, which are obtained by fitting 
the measured AI n ^\{V,<l>) and AI n: ± (V, <P) dependencies by Kulik's theory [65,66]. 



number n of Fiske step 


1/2 


1 


2 


3 


4 


5 


H\\ bottom electrode (x-direction): 


















140 


48 


27 


13 


6 


Qn\\ 




22 


27 


35 


30 


22 


Vn,\\ (MV) 


83.5 


167 


322 


475 


630 


772 


H _L bottom electrode (j/-direction): 














AI™? (jiA) 




130 


46 


24 


13 


5 


Qn,± 




19 


26 


31 


30 


18 


Vn,± QtV) 


90.5 


181 


336 


483 


636 


780 



junction versus the applied magnetic flux. For large n, the 
Fiske step has a maximum height at <P ~ n<P /2 where the 
Josephson current shows about the same spatial modula- 
tion along the junction as the cavity mode. 

The voltage positions V n of the Fiske resonances al- 
low us to derive several interesting junction parameters 
such as the Swihart velocity, the specific capacitance, or 
the junction plasma frequency. To avoid any ambiguities 
related to the overlap effect discussed above, we use the 
Fiske voltages found for H\\x to derive these parameters. 



Obviously, the values V n 



4>0£ 

2L' 



directly give the Swi- 



hart velocity c. From the position V\ = 165 /xV of the first 
Fiske step we obtain c to 0.027 • c, where c is the speed 
of light in vacuum. We note that the derived Swihart ve- 
locity may slightly vary with increasing step number, e.g. 
due to a nonvanishing dispersion of the dielectric constant 
of the barrier material [54,69]. Using the Josephson pen- 
etration depth Aj ~ 70 /im determined in section 3, the 
plasma frequency uj p /2it = c/Aj is obtained to 17.8 GHz. 
Since uj p = ^j2eJ c /hC s , we can derive the specific capaci- 
tance C s of the junction to C s = 54 fF /fim 2 . This value is 
about two times larger than typical values reported in lit- 
erature for Nb/AlOz/Nb Josephson junctions [70]. With 
these values the quality factor Qs W ihart = u} p R p AC s can 
be estimated, which directly follows from the Swihart ve- 
locity derived from the position of the first Fiske reso- 
nance. Here, R p is the resistance of the junction measured 
at a voltage V p — — 35 /LtV corresponding to the 

plasma frequency. R p agrees well with the subgap resis- 
tance i? sg = 0.44 fl. Using this value we obtain the qual- 
ity factor Qswihart = 6.6, which corresponds to the quality 
factor Qivc determined from the retrapping current of the 
IVC of the same junction using the RCSJ-model. The fact 
that Qswihart is slightly larger than the value Qivc =5.7 
is not astonishing keeping in mind effects of premature 
switching/retrapping due to noise in the IVC measure- 
ment. For the 20 x 20 /im 2 junction we observed Fiske 
resonances at V\ — 427 //V and V2 = 790 //V (not shown). 
Taking into account a slight increase of the first Fiske volt- 
age due to the overlap effect and a decrease of the second 
due to dispersion effects, we can estimate a Fiske step dis- 



tance of about 415 //V. This value leads to C s = 53 fF / /im 2 
in agreement with the value obtained for the large junc- 
tion. The derived plasma frequency is uj p /2it = 21.7 GHz 
and the quality factor is Qswihart =6.2. 

We next use the quality factors derived from the Fiske 
steps to analyze damping effects. In the theoretical de- 
scription of the Fiske resonances a finite damping is as- 
sumed. However, the origin of this damping is usually not 
specified. It turns out that the experimental values of Q n 
derived from the Fiske resonances can be both larger and 
smaller than the quality factor Qivc derived from the re- 
sistively and capacitively shunted junction (RCSJ) model. 
The reason is that on the one hand the RCSJ model over- 
estimates the losses due to quasiparticle tunneling, since a 
voltage independent resistance is assumed in this model. 
On the other hand, the RCSJ model does not take into 
account other loss mechanisms, e.g. due to a finite surface 
resistance. Assuming that there are various loss mecha- 
nisms the total quality factor of the Fiske resonances can 
be written as 



Qn — 



1 



1 



Qn 



qp 



Q, 



QthR s Qn,L 



\ 



i/Q„ 



(19) 



/ 



Here, the different contributions represent losses due to 
quasi-particlc tunneling (Q„ )q p) and a finite surface re- 
sistance (Q n ,R B ) as well as radiation losses (Q„, r ad), and 
dielectric losses (Q n ,e) [71]. Furthermore, variations AL 
in the junction length L lead to a broadening of the reso- 
nances which can be expressed by a frequency independent 
quality factor Q n ^ = I) /AD. 

If quasiparticle tunneling is the only damping mecha- 
nism we would expect Q n . qp = uj n RC. Here, R is the junc- 
tion normal resistance at the voltage V n = Huj n /2e, which 
is given by the about constant subgap resistance R sg in the 
relevant regime. Therefore, Q n .qp is expected to increase 
about linearly with increasing resonant mode frequency. 
Above we have determined Qswihart = w p i?C = 6.6. Ex- 
trapolating the value to larger frequency we would expect 
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Fig. 8. Dependence of the critical current I c and the height 
47™ ax of the Fiske steps on the magnetic flux generated by a 
magnetic field applied perpendicular (a) and parallel (b) to the 
bottom electrode. The data is obtained for a SIFS Josephson 
junction with A — 50 x 50 /im 2 at 500 mK. For clarity, the 
magnetic field dependence of I c and the height of the first 
Fiske step (Fl) are linked to the right axis. The data is fitted 
to a Fraunhofer pattern and Kulik's theory [65,66]. 



Qr 



0.37 • (w/27r)GHz . This dependence is shown 



in Fig. 9 by the broken straight line. Obviously, the qual- 
ity factors determined from the Fiske resonances do not 
follow this line. From this we can conclude that there are 
additional damping mechanisms beside quasiparticle tun- 
neling. In order to get some insight into the frequency 
dependence of these additional mechanisms, we have used 
(19) to determine the part in the quality factor due to the 
additional damping mechanisms by subtracting the quasi- 
particle damping. The resulting values Q n are shown as 
open symbols in Fig. 9. It is evident that Q n decreases 
significantly with increasing frequency. 

As discussed above there are several possible mecha- 
nisms leading to additional damping in Josephson junc- 
tions. However, radiation and dielectric losses usually can 
be neglected due to the large impedance mismatch at the 
junction boundaries and the small volume of the dielec- 
tric, respectively. Furthermore, in our fabrication process, 
I) /AD ~ 100 for I) = 50 resulting in Q n , L ~ 100. 
That is, the geometric inhomogeneities set an upper limit 



of the quality factor of the junctions, which is above the 
measured values. The remaining mechanism is damping 
due to the finite surface resistance of the junction elec- 
trodes. Expressing the complex surface impedance of a 
superconductor as Z s — R s + iX s , the contribution of the 
junction electrodes to the quality factor is given by the 
ratio [72] 



!n,R s = 



(20) 



For a rough estimate of Q n ,R s we can use a simple two- 
fluid model. Expressing the conductivity of the supercon- 
ductor as the sum a = a s + <r n of the superconducting and 
normal conducting carriers and using the simple relations 

[48] 



(21) 
(22) 



1 



- 10JT 



we can derive the following expression for the real and 
imaginary part of the complex surface impedance: 



Rs = 

2 



(^) 



(23) 
(24) 



Here, n n /n is the temperature dependent fraction of nor- 
mal electrons, m n the mass of the normal electrons, and 
<tn is the conductivity of the superconductor in the normal 
state where n n = n. With these expressions we obtain 



tn,R. 



u)/j, \1<j n (n D /n) 



(25) 



We see that Q n ,R B increases strongly with decreasing tem- 
perature due to the freeze out of the normal electrons. 
However, due to the uncertainties in n n /n, on in the 
PdNi/Nb bilayer, and the simplicity of our approach, 
eq. (25) certainly cannot be used to estimate the abso- 
lute value of Q n ,R a - However, according to (25) we expect 

Qn,R s oc l/u>. This explains the observed decrease of Q n 
with increasing frequency shown in Fig. 9. 

With the quasiparticle tunneling, the finite surface re- 
sistance and geometric inhomogeneities as the three main 
contributions to the measured quality factor, we expect 



1 



1 



qp 



Qn,R s Qn,L 



(26) 



As shown in Fig. 9, this expression well fits the measured 
data with Q° qp = 0.046 x 10" 9 s, Q° nJis = 2500 x 10 9 s"\ 
and Q n> L = 100. From this we can learn that the qual- 
ity factor of our SIFS junctions is limited by quasiparticle 
tunneling at low frequencies and the finite surface resis- 
tance at high frequencies. In the intermediate regime there 
may be an effect of geometric inhomogeneities when going 
to small area junctions. 



G. Wild et al.: Josephson Coupling and Fiske Dynamics in Ferromagnetic Tunnel Junctions 



13 



1 1 1 


1 ' 


i 


X < ' 1 ■ 

v> 


i 


- 


/o 






- 



















o 










♦ 
















€^ 








/'V 

/ / 














Q 

IVC 


■ 


1 / 




it 


Q 

Swihart 




1 / 






Q o Q 




.#' 




-■- 


Q i : Q 




/'.,., 






i.i. 


i 



50 100 150 200 250 300 350 



f (GHz) 

Fig. 9. Quality factors Qivc and Qswihart derived from the 
IVCs and the Swihart velocity by assuming only quasiparticle 
damping. The broken straight line gives an extrapolation of 
these values to higher frequencies. Also shown are the qual- 
ity factors Qii and Q± derived from the Fiske resonances (full 
symbols). The open symbols represent the quality factors Qii 
and Q± derived from the Fiske resonances after subtraction of 
the extrapolated quasiparticle damping. The dash-dotted line 
is obtained by fitting the data using (26). 



in the frequency regime between about 10 and 400 GHz 
achieved in our experiments are already sufficient for ap- 
plications in quantum information circuits or for studies 
of macroscopic quantum tunneling. 

We conclude the discussion of Fiske resonances by pay- 
ing attention to the small resonances labeled HF in Fig. 7, 
which occur at exactly half of the voltage of the first 
Fiske step. This observation reminds us of the appear- 
ance of half-integer Shapiro steps at the 0-7r-transition 
of SFS junctions due to a second harmonic component 
in the current-phase relation [76]. However, since the F 
layer thickness of the SIFS junctions studied in our work 
is about twice the thickness of the 0-7r-transition, this sce- 
nario would require a significant F layer thickness varia- 
tion of a few nanometers [77]. This is in contradiction to 
the small rms roughness of our F layers. Furthermore, the 
presence of a double-sinusoidal current-phase relation as 
the origin of the HF resonance is in contradiction with the 
measured magnetic flux dependence of its height. Clearly 
the measured oscillation period of the HF resonance can- 
not be mapped on twice the period of the first Fiske 
step [78]. Nonequilibrium effects may be a possible ex- 
planation for the observed HF resonance [79]. However, 
further experiments are required to clarify this point. 



5 Conclusion 



The quality factors measured for our SIFS Josephson 
junctions are slightly lower than the values reported for 
Nb/AlOz/Nb tunnel junctions [73]. This is not astonish- 
ing because there is the additional F layer in our SIFS 
junctions. First this layer reduces the I c R n product of the 
SIFS junction compared to a SIS junction and thereby in- 
creases the effect of quasiparticle damping characterized 
by Qivc - Qswihart = w p _R n C = ^2eI c RlC/h. Further- 
more, the F layer results in an increased surface resis- 
tance since the PdNi/Nb bilayer has increased ctn and 
an increased fraction n n /n of the normal electrons. The 
increased n n /n value is a result of the inverse proximity 
effect, increasing the quasi-particle density in the super- 
conducting electrodes [74] . Finally, effects originating from 
magnetic impurity scattering on Ni atoms diffused into the 
niobium top layer may play a role [75] . 

Comparing SIFS to SFS junctions, it is immediately 
evident that the quality factors of SFS junctions will be 
very small. Due to the much smaller normal resistance 
and vanishing capacitance, SFS junctions are overdamped 
(Q < 1). The larger quality factors of the SIFS junctions 
are obtained by the additional tunneling barrier, which 
causes large R n and C. Increasing the thickness t s of the 
tunneling barrier is expected to result in an exponential 
increase of R n , while I c R n should stay constant and C 
should decrease as Therefore, increasing t s in prin- 
ciple can be used to further increase the quality factors 
of SIFS junctions. However, this is obtained at the cost 
of lower J c and w p what may be a problem for some ap- 
plications. The quality factors between about 5 and 30 



We have fabricated Nb/Al(VPd . 82 Ni . 18 /Nb (SIFS) 
Josephson junctions with controllable and reproducible 
properties using a self-align process. High current den- 
sities up to more than 30 A/cm 2 and R n ■ A values above 
130 j?/im 2 have been achieved. The I c (^/<?o) dependencies 
are close to an ideal Fraunhofer diffraction pattern demon- 
strating the good spatial homogeneity of the junctions. 
The transition from 0- to 7r-coupled junctions was ob- 
served for a thickness — 6nm of the Pdo.s2Nio.1s layer. 
The Pdo.82Nio.1s layers show an out-of-plane anisotropy. 
They have a Curie temperature of 150 K, an exchange en- 
ergy E cx ~ 20 meV, and a saturation magnetization of 
about 1/xb per Ni atom, indicating that there are negligi- 
ble magnetic dead layers at the interfaces. The I c R n (d,F) 
dependence of the SIFS junctions can be well described by 
the dirty limit theory of Buzdin et al. [15] yielding the sin- 
gle characteristic length scale £f — 3.88 nm for the decay 
and oscillation of the critical current. From the measured 
IVCs and the Fiske resonances appearing at finite applied 
magnetic flux the junction quality factor has been deter- 
mined in the wide frequency range between about 10 and 
400 GHz. At low frequencies the quality factor increases 
about linearly with frequency due to the about frequency 
independent damping related to quasiparticle tunneling, 
whereas it decreases proportional to 1/w at high frequen- 
cies due to the increasing surface resistance of the junction 
electrodes. The achieved quality factors range between 
about 5 at the junction plasma frequency and 30 at about 
200 GHz and are sufficient for applications of SIFS junc- 
tions in superconducting quantum circuits or experiments 
on macroscopic quantum tunneling. 
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